Range of ion specific effects in the hydration of ions 



O 

C 
3 



Safir Merchant and Dilip AsthagirQ 

Department of Chemical and Biomolecular Engineering and The Institute of NanoBioTechnology, 
Johns Hopkins University, Baltimore, MD 21218 
(Dated: January 27, 2013) 

Within the quasichemical approach, the hydration free energy of an ion is decomposed into a 
chemical term accounting for ion specific ion-water interactions within the coordination sphere and 
nonspecific contributions accounting for packing (excluded volume) and long range interactions. The 
change in the chemical term with a change in the radius of the coordination sphere is the compressive 
force exerted by the bulk solvent medium on the surface of the coordination sphere. For the Na + , 
K + , F~, and Cl~ ions considered here this compressive force becomes equal for similarly charged 
ions for coordination radii of about 3.9 A, not much larger than a water molecule. These results 
show that ion specific effects are short ranged and arise primarily due to differences in the local 
ion-water interactions. 
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Numerous aqueous phase processes including, for ex- 
ample, the stability of biological macromolecules and col- 
loidal suspensions, are strongly influenced by the chem- 
ical type of the inorganic ions dissolved in the aqueous 
phased. However, despite their ubiquity, a molecular level 
understanding of ion specific effects remains elusive^. 
Ideas such as distortion of hydrogen bond network of 
the bulk solvent^., balance between ion-water and water- 
water interactions^, and influence of ion-specific disper- 
sion forces^ have been proposed to rationalize these ef- 
fects. 

Understanding how the ion affects the solvent ma- 
trix remains of first interest in understanding ion spe- 
cific effects. Recent spectroscopic measurements on salt 
solutions^- 8 - suggest that the ion only affects hydro- 
gen bonds of water molecules within its first hydration 
shell. An earlier theory and simulation study by us^ also 
showed that only a small subset of water molecules in 
the first hydration shell of the ion are sensitive to the 
type of the ion. These results taken together suggest 
that ion specific ion-solvent interactions are limited only 
to the ion's local neighborhood. In this Communication, 
we address the question of how far in the liquid the chem- 
ical type of the ion is felt and if the point at which ion- 
specificity is lost is also the point where continuum mod- 
els of hydration begin to take hold. 

To address these questions, we first define a coordina- 
tion sphere of radius r around the ion to separate the ion 
interaction with water molecules within the coordination 
shell from the longer-range interaction of the ion with 
the solvent outside the coordination sphere. With such a 
spatial partitioning, the hydration free energy of the ion, 



/i cx , can be written asi£ — 

H cx = k B T In x (r) - k B T \n p Q (r) 



Pouter ( r ) 



(1) 



where fee is Boltzmann's constant and T is the temper- 
ature, k B T 'In xo(r) is the chemical contribution to the 
hydration free energy, — k B T\ripQ(r) is the packing con- 
tribution, and Pouter ( r ) i s the contribution to the hydra- 
tion free energy due to ion interaction with the solvent 
outside the coordination shell. Xo(r) is the probability 



of observing no water molecules within the coordination 
sphere of radius r around the ion; thus k B T\rixa{r) is 
the free energy gained by allowing water molecules to 
flood the empty coordination sphere, justifying the iden- 
tification of this term with local, chemical interactions. 
po (r) is the probability of observing an empty coordina- 
tion sphere, but in the absence of the ion; — k B Thipo(r) 
is the free energy required to create a cavity of radius r 
in the bulk solvent. Mouter( r ) ^ s the free energy gained 
in placing the ion in an empty cavity; it accounts for the 
interaction of the ion with the bulk medium outside the 
coordination sphere. For a sufficiently large r, we ex- 
pect Pouter ( r ) to be composed of a large number of small 
non-specific contributions^^. Since Eq. [1] is a tautol- 
ogy, the point where MouterM becomes insensitive to ion 
type is also the range to which /cbTIii^oM captures all 
the ion-specific effects. 

To identify the range of ion specific effects, we studied 
ion- water systems under NVT conditions using Metropo- 
lis Monte Carlo simulation o 16 i 17 . Water was modeled us- 
ing the SPC/E model 18 and parameters for Na + , K + , 
F _ , CI - ions were obtained from an earlier studyi^. The 
chemical term was estimated by growing the coordination 
sphere in steps of 0.05 A. At any r, the probability of ob- 
serving zero (0) water molecules in a coordination sphere 
of radius r + 0.05 A was obtained. From this data, the 
entire xq(t) versus r profile was reconstructed. A simi- 
lar strategy was followed for po(r) versus r, but with a 
step size of 0.1 A. At each r value, 3 x 10 5 sweeps of 
equilibration were followed by 3 x 10 5 sweeps of produc- 
tion. Each sweep comprised one attempted move of each 
particle in the system. The maximum allowed transla- 
tion and rotation were adjusted during the equilibration 
phase to target an acceptance ratio of 0.3. In the pro- 
duction phase these maximum displacement values were 
held fixed. Configurations are saved every 10 sweeps for 
analysis. 

From Fig. Q] we observe that for r > 3.9 A the change 
in the chemical term depends only on the charge of the 
ion and is independent of its chemical nature. Fig. Q] 
emphasizes that, for the ions studied here, ion-specific 
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FIG. 1: The change in the chemical term, AksTlnxo, with 
change in the radius, r, of the coordination shell for various 
ions. The dashed line is the result for a solute that does not 
interact with the solvent; it thus accounts for the change in 
the packing contribution. Observe that the change in the 
chemical term approaches the change in the packing term as 
the coordination shell is increased. For r > 3.9 A the change 
in the chemical term is the same between identically charged 
ions. 

effects are short ranged. In particular, for r > 3.9 A, 
the increase in the chemical term is due to non-specific 
contributions. 

To elucidate the non-specific nature of contributions 
for r > 3.9 A, we first note that 

d\nxo(r)/p (r) d^Z tcr {r) 
kBT dr = — ■ (2) 

The left hand side of the above equation is the compres- 
sive force of the bulk on the coordination sphere of the 
ion relative to the compressive force on an ion-free coor- 
dination spher 

For a sufficiently large r, we expect the interaction of 
the ion with the bulk medium outside the coordination 
sphere to be well-described by a Gaussia n 14 i 15 . In this 
case, for an ion of charge q, we have 

= ?Wo-f -g 2 -((0-(0)o) 2 }o, (3) 

where (. . .)o denotes averaging in the presence of an 
uncharged ion, {cf>)o is the potential at the center of 
an uncharged ion, and {{(f) — {4>)o) 2 )o is the fluctuation 
in the electrostatic potential^. By defining 1/rBom = 
(3{(4> — (0)o) 2 )o/2, it is readily seen that the fluctuation 
contribution is just the Born model of hydration. (For 
clarity, in Eq. [3] we do not explicitly indicate the cor- 
rections due to electrostatic self-interaction-^, as these 
depend solely on q and not on the ion type.) 

In the radius range of interest we find that d{4>)o/dr 
is about a k-gT / 'A. (A similar estimate is obtained based 
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FIG. 2: F = k B T aiux ° { g ) r /po{r) , Eq. H the compressive force 
exerted by the solvent on the surface of the coordination shell 
containing the ion, k^Td\\n xo{r)]/dr, relative to the com- 
pressive force on an empty cavity, kBTd[\npo{r)]/dr for var- 
ious coordination radii, r. The black dashed line is the es- 
timate of the relative compressive force based on the Born 
model of hydration (Eq.2}. 

on the data in Ref. 21 for a cavity.) Neglecting this con- 
tribution, we thus expect that asymptotically 

dlnxo/p q 2 1 
kBT — 2^ {1 --e ] > (4) 

where e is the dielectric constant for water. 

Figure [2] shows that for r > 3.9 A the relative compres- 
sive force between similarly charged ions becomes equal, 
but is different for anions and cations. The large r value 
is also different from the value obtained using the Born 
model, where the Born-radius is set equal to the coordi- 
nation radius. 

The discrepancy in the asymptotic limit between sim- 
ulation data and that based on the Born model is not 
particularly surprising, since the coordination radius is 
not the same as the Born radiu s 19 i 21 . For example, the 
position of the water molecule is denoted by the position 
of its oxygen atom and thus the coordination shell strictly 
excludes only the oxygen atoms of the water molecules 
but not their corresponding hydrogen atoms, which can 
enter the coordination sphere. Since this effect will be 
more pronounced for the anions than the cations, we ex- 
pect that the equivalent Born radius for anions to be 
smaller than that for the cations. What is interesting, 
however, is the observation that by empirically reduc- 
ing the cavity radius by 1.3 A for anions and 0.5 A for 
cations (Fig. [3]) , the relative compressive force for both 
anions and cations are similar for large coordination radii. 
Further, the large-r curvature of the radially translated 
curves compares reasonably well with that based on the 
the Born model of hydration. (We maintain the differ- 
ence of 0.8 A, as earlier studies on such shifts find a sim- 
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FIG. 3: The data in Fig.[2]have been shifted along r by —1.3 
A for the anions and by —0.5 A for the cations. With this 
empirical shift, the compressive force F — ks r ^n^oM/poM ^ 
Eq. [21 follows the trend predicted by the Born model (Eq. [4]| . 



In summary, for the ions studied here, we find that a 
coordination radius of about 3.9 A, not much larger than 
the size of a water molecule, is sufficient to account for all 
ion specific ion-water interactions. Our results are in ac- 
cordance with experimental results^— that suggest that 
the influence of the ion on hydrogen bonding between wa- 
ter molecules is minimal for water molecules beyond the 
first hydration shell. The short-range of ion specific ion- 
water interactions further suggests that any framework 
for modeling ion-specific effects needs to acknowledge the 
molecular characteristics of local ion-water interactions. 
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